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Abstract. The recent X-ray and CO observations of RX J1713.7—3946 show that a significant 
fraction of the nonthermal X-ray emission of this unique supernova remnant associates, in one way 
or another, with a molecular cloud interacting with the west part of the shell. This adds a new 
puzzle in the origin of X-ray emission which cannot be easily explained within the standard model 
in accordance of which X-rays are result of synchrotron radiation of multi-TeV electrons accelerated 
by strong shock waves. We explore an alternative origin of the X-ray emission assuming that it is 
produced by secondary e* resulting from high energy hadronic interactions in the molecular gas. 
Such a scenario could explain in a quite natural way the apparent correlation between the X-ray 
and CO morphologies. However, the TeV y-ray emission recently reported by H.E.S.S. significantly 
constrains the parameter space of this model. Namely, this mechanism cannot reproduce the bulk of 
the observed X-ray flux unless one postulates existence of a PeV cosmic-ray component penetrating 
with an unusually hard spectrum into the dense cloud. 


SNR RX J1713.7—3946: BRIEF SUMMARY 

The SNR RX J1713.7—3946 exhibits remarkable nonthermal X-ray emission, with the 
largest total energy flux of ~ 5 x 10~ 10 erg cm~ 2 s -1 (2-10 keV) amongst the shell- 
type SNRs. The diffuse X-ray emission, presumably of synchrotron origin, is distributed 
over the entire SNR. It is characterized by a power-law spectrum with an average photon 
index 2.3 

The observations of this source with Chandra 0,1 revealed that the bright north¬ 
western (NW) rim is comprised of a complex network of filaments and hotspots em¬ 
bedded in diffuse emission. The X-ray spectra of these sub-components in the NW rim 
are well fitted by a power law with variation of the photon indices within T ~ 2.1-2.5 
OD - Recently, similar spatial and spectral characteristics have been found in the south¬ 
western (SW) rim in the XMM-Newton data 0i. The applicability of the standard 
shock-acceleration (more specifically, diffusive shock acceleration) scenario to the syn¬ 
chrotron X-ray emission in RX J1713.7—3946 has been disputed based on the lack of 
indication of the expected (within the framework of this model) cutoff energy below 10 


keV|H. 

Originally it has been claimed that the molecular clouds found at a distance of 6 
kpc are associated with this SNR |2]. Therefore the distance to the SNR has been widely 
adopted as 6 kpc. However, the recent CO observations made with the sensitive NANTEN 
telescope have revealed a strong indication of interaction of a molecular cloud located at 
a distance of 1 kpc with RX J1713.7—3946 0]. This finding is further strengthened by 
an elaborate comparison between the X-ray and molecular emissions presented below 
and by the X-ray absorption features [8]. The atomic and molecular hydrogen column 
density toward the line-of-sight of the X-ray nebula becomes comparable to X-ray 
absorbing column density at ~ 1 kpc [9], which also supports the association with the 
nearby molecular cloud. The ~ 1 kpc distance to the supernova remnant makes plausible 
its association with the guest star A.D. 393 [ 10] (the remnant age of / age — 1600 yr). 

The CANGAROO collaboration reported the detection of TeV y-ray emission from 
the direction of the NW rim EUH. The spectrum was claimed to be quite steep with 
a power-law photon index F = 2.8 ±0.2 (0.4-8 TeV). It was argued that the steep 
spectrum is inconsistent with the inverse Compton model, but could be explained by 
7T°-decay y-rays 0. Most recently, the H.E.S.S. collaboration announced detection of 
TeV y-rays distributed over the entire X-ray nebula with a hard power-law spectrum 
r = 2.2 ±0.1 (1-10 TeV) and flux 3.5 x 10“ 11 erg cm~ 2 s -1 iHcll . The fascinating y-ray 
image detected by H.E.S.S. appears to be generally similar to (although not exactly same 
as) the X-ray image. 

In this paper we compare the spatial structures of the nonthermal X-ray nebula 
RX J1713.7—3946 with the CO image of the interacting molecular cloud, and briefly 
discuss the origin of X-radiation which seems to be quite different from synchrotron 
X-ray emissions of other young shell-type SNRs (14], and explore the capability of 
secondary e ± , generated through hadronic collisions in the cloud, to explain the X-ray 
nebula. Details will be presented elsewhere [8|. 


RX J1713.7—3946: COMPARISON BETWEEN NONTHERMAL 
X-RAY AND CO EMISSIONS 

In Fig. |H the overall X-ray nebula is compared with the NANTEN CO distribution of 
the nearby interacting molecular cloud. The X-ray nebula shapes a large parallelogram 
(60' x 70 / ), with enhanced brightness in the northwest and southwest parts, which we 
refer as “NW rim” and “SW rim”, respectively. These rims with an apparent width 
~ Kf have been extensively examined in Refs. [3:1 HI EL 6]. The bright NW and SW 
rims appear to be buried in the molecular cloud. Four CO peaks (A, B, C, and D) are all 
located in the western part of the nebula. A direct correlation between the X-ray and CO 
emissions can be seen for peak C, from which the broad-line component of CO emission, 
a good indicator of SNR interactions 0 , has been detected |@]. This morphological 
comparison suggests that the dense molecular gas somehow plays a crucial role for the 
formation of the bright nonthermal X-ray rims. 

A close-up X-ray image of the SW portion, where the major CO peaks are located, 
is displayed in Fig. |2ji together with the CO(7=1-0) intensity map. The XMM-Newton 
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FIGURE 1. ASCA X-ray image (1-5 keV) of RX J1713.7 3946 overlaid with CO(J=1-0) map 
(contours) in the velocity range from —11 to —3 km s -1 . The X-ray intensity scale is counts ks 1 pixel -1 
with a pixel size of ~ 0.5'. The CO contours are every 2 K km s 1 in the range of 2-30 K km s 1 . North 
is up and east is to the left. 



FIGURE 2. XMM MOS images (0.6-5 keV) accumulated in 10" x 10" pixels with the same CO map 
used in Fig.Qoverlaid. The X-ray intensity scale is counts ks -1 pixel -1 . North is up and east is to the left, 
a) For the SW part. The CO contours are every 2 K km s -1 from 2 to 30 K km s. b) For the SE part. 
The CO contours are from 1 to 8 K km s -1 with an interval of 1 K km s. 

















image demonstrates that the CO peak C has excellent spatial match with the nonthermal 
X-ray emission. Through spectral analysis, we found that the X-ray emission from peak 
C experiences additional soft X-ray absorption (A/Vh — 2 x 10 21 cm -2 ), confirming the 
association with the CO peak. Also we found another bright X-ray feature that correlates 
well with the CO clumps, lying at the northeastern edge (a steep gradient of the CO line 
intensity) of peak A (A1+A2). The X-ray spectra of these bright features have a photon 
index T ~ 2.3, quite typical for the west part of the rim. 

In Fig. |2t>, we also present the X-ray image of the southeastern (SE) sub-field, where 
the X-ray intensity is very low. The faint diffuse X-ray emission seems to associate 
with the molecular gas in the SE field. Interestingly, while in the SW rim the brightest 
X-ray features emanate from the molecular cloud, in the SE part the faintest appear to 
originate from the cloud. The X-ray spectrum extending up to at least 5 keV suggests 
nonthermal origin, although the faintness of the X-ray emission from this region does 
not unfortunately allow a detailed spectral study. 

The results presented in Figs. Sub imply that the nonthermal X-ray emission emanates 
directly from dense molecular gas (though it cannot be excluded that a part of the emis¬ 
sion originates in hot tenuous plasma). Usual synchrotron interpretation for the X-ray 
emission requires presense of multi-TeV electrons in the molecular cloud. Acceleration 
of multi-TeV electrons inside the molecular cloud, however, is unlikely with conven¬ 
tional shock-acceleration models. Slow shock waves propagating inside the dense cloud 
(v <C 1000 km s _1 ) are not fast enough to accelerate electrons to very high energies; the 
diffusive shock acceleration time scales as f acc v~ 2 , where v is the shock velocity. 

Yet external acceleration scenarios also involve some difficulties. Suppose, for ex¬ 
ample, that the fast supernova shock expanding into the stellar wind bubble accelerates 
multi-TeV electrons and injects them into the molecular cloud resulting in synchrotron 
X-ray emission. The ultra-relativistic electrons are subject of severe synchrotron losses 
as passing through the magnetized cloud, with synchrotron cooling time f syn ~ 240 yr 
for a moderate magnetic field of B — 30 /lG. Unless the X-ray nebula is extremely young 
(f a ge < fsyn), one should observe aging effects of synchrotron-emitting electrons. Never¬ 
theless, we measured hard synchrotron X-ray spectra (r ~ 2.3) throughout the NW and 
SW rims, which seems to argue against external origins of the X-ray emitting electrons. 

The power-law spectra of synchrotron X-rays measured in this nebula suggest syn¬ 
chrotron cutoffs located at energies beyond £ 0 10 keV |@J. On the other hand, within 

the standard framework of diffusive shock acceleration, the synchrotron cutoff should 
be located well below 10 keV. 

Therefore it is important to explore other possible (alternative) models for interpre¬ 
tation of the spectral and spatial features of the observed nonthermal X-ray emission. 
Below we discuss the capability of a hadronic model to explain a non-negligible fraction 
of X-rays by synchrotron emission of electrons from decays of secondary 7t ± mesons 
produced in hadronic interactions. 


SECONDARY-ELECTRON SYNCHROTRON X-RADIATION 


Interactions of accelerated protons with the ambient gas lead to the production of 
7T-mcsons with comparable fractions of the energy of primary protons transferred to 
the final products of /r-meson decays—y-rays, electrons and neutrinos. While y-rays 
and neutrinos freely escape their production regions and can be observed directly, 
the secondary electrons lose their energy through synchrotron radiation and inverse 
Compton scattering. In presence of sufficiently large magnetic field, e.g. close to 100 /iG 
or more, which is needed for acceleration of protons to PeV energies, the bulk of the 
electron energy is released in the form of hard X-rays. The production of X-rays by 
secondary electrons is an attractive possibility which to a large extent is free of the 
problems associated with the synchrotron radiation by directly accelerated electrons. 
In particular, this model allows diffuse synchrotron X-ray emission throughout the 
remnant (because the protons can propagate without significant losses), as well as allows 
extension of X-ray spectra well above 10 keV. 

One may expect that directly accelerated electrons are dominant contributors to syn¬ 
chrotron X-ray emission, rather than the secondary electrons. Actually, indeed the pri¬ 
mary synchrotron X-radiation can exceed the contribution by the secondary electrons 
by more than an order of magnitude for shell-type SNRs with low ambient density, e.g. 
n ~ 1 cm~ 3 and a proton-to-electron ratio ~ 100. However, since the cease of particle 
injection quickly results in a dramatic drop of the X-ray flux by primary electrons on a 
timescale / syn ~ 100 yr, while the X-ray flux by secondaries remains unchanged on much 
longer timescale T pp <; 10 5 yr, in the case of a “dead” accelerator (e.g. due to the crash of 
the shell with a dense cloud) we may expect X-radiation only from secondary electrons. 
Moreover, in dense environments, like molecular clouds interacting with the SNR shell, 
the synchrotron X-radiation may dominate even in the case of an active accelerator. 

Finally, the long-lived ultrarelativistic protons can propagate throughout the dense 
gas generating secondary e ± via the production and subsequent decay of charged pions. 
Unlike directly accelerated electrons which are expected to be concentrated in the 
proximity of their acceleration sites, the secondary electrons can be easily distributed 
over larger volume, resulting in extended synchrotron radiation as observed. 

In Fig.|3] we show a characteristic example of broad-band SEDs of nonthermal radi¬ 
ation resulting from pp interactions, assuming that protons are injected into the produc¬ 
tion region with a (number) spectrum Q °e E~ l 9 e~ E ^ E ° . It is seen that for the sufficiently 
high cutoff energy in the spectrum of protons (Eq > 1 PeV), the secondary synchrotron 
X-ray flux becomes comparable to that of 7T°-decay y-rays. Since the 7T°-dccay y-rays 
are tightly coupled in such a way with the synchrotron radiation of secondary electrons, 
the TeV y-rays flux detected by H.E.S.S. [ 13] constrains the level of the X-ray flux 
contributed by secondary electrons. Note that the n° components shown in Fig. [3] cor¬ 
respond to the H.E.S.S. flux from the NW part. Since the accompanying X-ray flux by 
secondary-electrons falls well below (by an order of magnitude) the observed X-ray flux 
(at several mCrab level), one may conclude that the contribution by secondary electrons 
cannot exceed 10 percent of the observed X-ray flux. 

However, the relative contribution of X-rays can be significantly increased without 
violating the H.E.S.S. flux if one assumes very hard proton spectrum, e.g. Q °c E -1 - 5 , 
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FIGURE 3. The broadband SEDs of radiation (secondary e ± synchrotron and 7T°-decay gamma-rays) 
initiated by p-p interactions. It is asumed that protons are injected continuously into a dense region with 
linear size R = 1 pc, gas density n = 100 cm 3 , and magnetic field B = 100 uG. The injection spectrum 
of protons is a power law with spectral index s = 1.9 and exponential cutoff Eq = 0.1 PeV, 1 PeV, and 
10 PeV. The source age 1000 yr. The proton injection power L p = 3 x 10 37 erg s 1 (total injected energy 
W p ~ 1 x 10 48 erg). The escape time of protons in the source is set assuming Bohm diffusion. 


extending to PeV energies. It should be noted that such a hard component of cosmic rays 
is not excluded by the nonlinear shock acceleration theory ES. Also, very hard proton 
spectra can be expected due to propagation effect. For example, it is possible, at least in 
principle, that due to very slow diffusion only > 1 PeV protons can arrive, for the limited 
age of the accelerator, at the dense regions of the clouds. In this case formally we can 
expect an arbitrary hard proton spectrum which in the regions of enhanced density may 
initiate secondary radiation with unusual SEDs—with distinct peaks at the hard X-ray 
and ultrahigh energy y-ray domains. 
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